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Abstract 

The perovskite nickelates RN 1 O 3 (R: rare-earth) have been studied as potential multiferroic compounds. A cer¬ 
tain degree of charge disproportionation in the Ni ions has been confirmed by high resolution synchrotron power 
diffraction: instead of the nominal Ni 3+ valence, they can have the mixed-valence state Ni (3_3)+ and Ni (3+<>)+ , though 
agreement has not been reached on the precise value of 6 (e.g. for NdNi 03 , 6 = 0.0 and <5 = 0.29 were reported). 
Also, the magnetic ground state is not yet clear: collinear and non-collinear Ni-O magnetic structures have been 
proposed to explain neutron diffraction and soft X-ray resonant sccattering results in these compounds, and more 
recently a canted antiferromagnetic spin arrangement was proposed on the basis of magnetic susceptibility mea¬ 
surements. This scenario is reminiscent of the situation in the half-doped manganites. 

In order to gain insight into the ground state of these compounds, we studied the magnetic excitations of some 
of the different phases proposed, using a localized spin model. With the purpose of describing the charge dis¬ 
proportionation, we include two kinds of Ni-spins with different magnitudes. As for the magnetic couplings, we 
include: nearest-neighbor (NN) and next-nearest-neighbor (NNN) Heisenberg-like interactions, respectively for 
the ferromagnetic and antiferromagnetic couplings present in the collinear phases. To describe the non-collinear 
phases, and as already proposed for other multiferroics, we also consider NN Dzyaloshinskii-Moriya-type cou¬ 
plings to allow for the possibility of a relative angle 6, between NN spins in the two different magnetic sublattices. 
Using a simplified spin chain model for these compounds, we first analize the stability of the collinear, orthogonal, 
and intermediate phases in the classical case. We then explore the quantum ground state indirectly, calculating 
the spin excitations obtained for each phase, using the Holstein-Primakoff transformation and the linear spin-wave 
approximation. For the collinear and orthogonal (6 = n/2) phases, we predict differences in the magnon spectrum 
which would allow to distinguish between them in future inelastic neutron scattering experiments. 
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1. Introduction 

The ferroelectric oxides with magnetic ordering have attracted much attention since they offer the possibility of 
controlling the electric polarization or the magnetic ordering by applying magnetic or electric fields, respectively, a 
desirable feature in the design of electronic devices (Ij]. However, finding these multiferroic oxides has not been an 
easy task. Though there are some of them which have a simultaneous ferroelectric character and magnetic ordering, 
usually the coupling between these is very weak and therefore poorly controlled with applied fields. In 2004 
Efremov et al. yfl suggested that in manganites (RMnC >3 R: rare-earth), in addition to simultaneous charge and 
magnetic ordering, a charge disproportionation (CD) of the Mn ions would be needed for these materials to become 
multiferroic. This CD means that instead of the nominal valence Mn 3+ , mixed valences Mn l3 ~ li)+ and Mn (3+i)+ 
should be present. Later, van den Brink and Khomskii JjJ discussed about the possibility of ferroelectricity related 
to charge disproportionation in rare earth perovskite nickelates of the type RNiC >3 (R=rare earth). In fact, in 2000 
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Mizokawa et al. 0] had studied a multiband d - p model for perovskite transition metal oxides, suggesting that 
it could describe those nickelate compounds, and found an antiferromagnetic ground state with charge ordering 
centered either in the 0-2p orbitals, for relatively large charge-transfer energy (as in PrNi 03 and NdNiCT), or with 
charge-ordering in the transition metal 3d orbitals, for negative charge transfer energies (relevant for YNiOs). 

These nickelates (7?=rare-earth, or Y) present a metal-insulator transition 00 ] at temperature Tmi, and anti¬ 
ferromagnetic ordering below the Neel temperature 7\(^ Tmi) with a possible ordering of Ni (3- ' 5)+ and Ni (3+l5)+ 
ions i8-171. with various values of the charge disproportionation 6 as discussed below. As one example, in NdNiC >3 
it was found that Tmi = 7’n = 200 K 0 , 0 . 

In 2009, Giovanetti et al II18fl showed by first principles calculations, that in nickelates simultaneous charge 
and magnetic ordering could be present, as well as a charge disproportionation of the Ni ions, and electrical 
polarization would thus be induced. In their work, they calculated the electrical polarizations obtained for three of 
the magnetic phases previously proposed for nickelates, along with a specific charge ordering of Ni 2+ and Ni 4+ ions 
corresponding to a charge disproportionation of 5 = 1. The magnitude and direction of the electrical polarization 
induced would indicate the underlying magnetic order in these oxides, at present still not clear. 

As shown in Figure [T| the magnetic orderings which they studied 1 11811 are: i) the S-collinear phase first 
proposed by Garcia et al. 18[] in X-ray (XRD) and neutron diffraction (ND) experiments for PrNiOs and NdNi 03 , 
and later by Fernandez et al. JB for HoNi 03 ; ii) the T-collinearphase proposed by Giovannetti ifisll : and the iii) 
N-non-collinearphase proposed by Scagnoli et al. EG!] for NdNi 03 based on soft X-ray resonant scattering at 
the Ni-Lo.s and Nd-M edges. 
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Figure 1: (Reproduced from Ref Ml 8 ll -Fig.l.) Schematic view of the charge and magnetic structures of RMO 3 
oxides, (a) Collinear up-up-down-down magnetic structure 0], (b) Non collinear magnetic structure 1131 


The S-collinear phase (see Figure |T). is characterized |fl~ 8 ll by a checkerboard charge order of Ni (3 ~ l5)+ and 
Ni (3 l0) 1 ions, corresponding to spins S 1 and S 2 respectively, along with a magnetic structure defined by the propa¬ 
gation vector k = (1/2,0,1/2), not seen in other perovskite oxides. This involves alternating ferromagnetic (FM) 
and antiferromagnetic (AF) couplings along the three pseudocubic axes such that every Ni-spin is coupled FM with 
three of its nearest neighbors (NN) and AF with the remaining ones. Regarding the magnetic cell, this structure 
can be pictured as formed by ab planes stacked in c direction in the form A 1 A 1 A 'A", where in A~ all spins are 
inverted with respect to A + . Notice that on each plane, there are FM zigzag chains along b, which are coupled AF to 
each other. Experimentally, the direction of the moments within each plane appears to be either along a (Ref. 10), 
or in the ac plane (Refs. |JU). J_2|). Notice that the T-collinear phase , differs from the S-phase in the stacking of 
the zig-zag chains between adjacent planes: in the S phase all zig-zag chains point in the same direction, whereas 
in the T phase in alternate planes they point in o ppo site directions Il 8 ll . The N-non-collinear phase has the same 
charge order as both collinear phases considered |[l 8 |]. However, its magnetic structure corresponds to a spin spiral, 
in which the spins in FM planes perpendicular to the [101] direction appear rotated around the [010]-axis between 
consecutive planes. Note that this N-non-collinear phase is different from other non-collinear phases proposed for 
nickelates: in Ref. 112], these planes are alternatively FM and AF, while in Ref. E the FM ab planes are stacked 
and rotated 6 ~ 76° along [001]. Apart from these phases, recently a canted antiferromagnetic spin arrangement 
was suggested on the basis of magnetic susceptibility measurements [ 20 h. 

A wide set of values has been reported for the charge disproportionation 5 found in different rare-earth nicke¬ 
lates, as we describe next. For PrNi 03 and NdNi 03 , the first studied compounds, 5 = 0 according to Refs. 00 . 
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However, more recently 6 ~ 0.21 was reported for PrNiC >3 [14|], whereas for NdNiC >3 a value of 6 ~ 0.29 was 
estimated in Ref. [16] while Ni l ' 2 - 5±,i ’ l+ states with o' ~ 0.16 follow from Ref. [13]. For YNi 03 in Ref llToh 6 ~ 0.28 
was estimated, which coincides with the value in the study through the whole series of R =Y, Ho, Er, Tm, Yb, Lu 
in Ref. liTHl . where 5 =0.28, 0.38, 0.32, 0.36, 0.33, 0.33, were respectively reported. For HoNiC> 3 , nevertheless, a 
larger value 5 ~ 0.48 can be estimated from the reported magnetic moments in Ref Q. ForTmNi 03 and YbNiC> 3 , 
from isomeric shifts in Ref. 5 ~ 0.14, and 0.16 respectively, were estimated, values which correspond to ap¬ 
proximately half the indicated CD in Ref tm. For DyNiC >3 in Ref. [G~§] 6 ~ 0.52 is found for the non-collinear 
phase which best agrees with their experiments. 

In the present work, as a first approach to the study of the problem in nickelates, we study the magnetic 
excitations of a one-dimensional (ID) chain, like the ones included in the collinear and non-collinear phases 
analized in Ref. Q. We used the localized spin model to be presented in next section, where the possibility of 
charge disproportionation is included by considering Ni-spins with eventually different magnitudes. Regarding 
magnetic couplings, in our model we include the minimal set required to describe the collinear as well as the non- 
collinear phases proposed. That is, FM nearest-neighbor (NN) and AF next-nearest-neighbor (NNN) Heisenberg 


like interactions, and to describe non-collinear phases also a NN Dzyaloshinskii-Moriya-type (DM) coupling [21 


12211 . to allow for the possibility of a relative angle 6 between NN spins in different magnetic sublattices. We analize 
the stability of the collinear, orthogonal, and intermediate phases in the classical model. Then, we explore the 
quantum ground state indirectly, by calculating the spin excitations obtained for each phase, showing that for the 
collinear and orthogonal (6 - n/2) phases proposed, differences in the magnon spectrums are to be expected, 
which would allow to distinguish between them experimentally. 


2. Simplified chain model: generic intermediate phase 

In order to describe the main ground state phases proposed for nickelates mentioned in previous section, we 
study a simplified chain model and propose a generic phase, which we call the “intermediate phase” which, as 
respective limiting cases, can describe the collinear and the orthogonal phases. As a first approach, here the three- 
dimensional nickelate compounds are studied using a simplified model: representing them by chains of localized 
spins, shown by dashed lines in Figure l2al 



(a) (b) 

Figure 2: (a) Schematic view of the charge and magnetic orderings under discussion for RNiC >3 oxides. With 
dashed lines we represent the spin chains present, object of our present study. The charge disproportionation 
Ni (3 0)1 and Ni (3+<5,+ is represented, respectively, by gray and black circles, (b) Generic intermediate phase, for one 
spin chain: i.e. a spin wave, with spins of alternate magnitude. Here the S i-sublattice is considered fixed, while 
the S 2 -sublatttice is rotated by an angle 6 with respect to the former. 


In the localized spin model we consider two kinds of Ni-spins with different magnitudes in order to describe the 
charge disproportionation. As shown in Figure[2b] and to take into account the phases proposed in experiments 
10, 12, 3, along the chains we consider a unit cell composed by four spins: two of them with magnitude S \ 


representing the Ni (3 3,+ sublattice, and two other ones with magnitude Si for Ni (3l0)l ', being 6 a measure of 
the Ni charge disproportionation (CD). The main difference between the experimentally proposed collinear and 
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orthogonal phases is the relative orientation between the two antiferromagnetic sublattices, which we describe by 
angle 0 characterizing the intermediate phase, shown in Fi gure [Zbl Notice that the sign of 6 will determine the 
“helicity of the spin chain. The collinear phase jg, kJ 12, Tstl is characterized by 9 = 0, while the orthogonal 
phase | jjl TU corresponds to 9 - zr/2 and all chains have equal helicity. It is worth mentioning that in Ref. 1121a 
slightly different non-collinear phase is proposed: with 6 ~ 0.447T, and anisotropic helicity (along the z-direction, 
alternating “helicity is proposed for consecutive chains). 

Regarding magnetic couplings, in our model we include: FM nearest-neighbor (NN) and AF next-nearest- 
neighbor (NNN) Heisenberg-like interactions, to describe the collinear phases. To describe non-collinear phases, 
the model also includes a NN Dzyaloshinskii-Moriya-type (DM) coupling |2ll 22], like the one previously used in 
Ref. 112311 to analize the spin excitations in the distorted NiOi planes in La 2 Ni 04 , in order to allow for the possibility 
of a relative angle 9 between NN spins in the two different magnetic sublattices. We found that the minimal model 
of localized spins which could describe the single chains present in the two limiting phases experimentally pro¬ 
posed HOI as well as the generic intermediate phase with other 9 values as in Ref. 112}] , requires the inclusion of 
those three magnetic couplings. In particular, if we picture the 4-spin unit cell as formed by two “dimers” (plotted 
in different color/linewidth in Figure [2bl>. an “intra-dimer” ferromagnetic (FM) coupling F and a Dzyaloshinskii- 
Moriya (DM) like coupling K are required, as well as a NNN antiferromagnetic (AF) coupling A in each magnetic 
sublattice (see FigurelZbl). 

With the above considerations, we studied the following spin chain Hamiltonian: 


<H = 


-F 2, V 

( n,m)/eD 


z 


+ K > y 

(n,m)/eD 


(S n x S'j + A ^ S« 

in,ml) 


( 1 ) 


where D indicates spins inside dimers, (n, m) or (( n,m )) indicate nearest-neighbor (NN) or next-nearest-neighbor 
(NNN) spins, respectively. Here, all couplings are considered positive, and the sign of K determines the helicity 
of the spin chain: K > 0 produces a counterclockwise rotation like in Figure [2b] while K < 0 would induce a 
clockwise rotation. The primes refer to the use of 0-rotated local spin quantization axes for S 2 -spins, as detailed 
later. 


3. Results and discussion. 

3.1. Phase stability in the classical model 


As a first step, we analize the classical behaviour of the spin chain model of Equation £[]). The energy of the generic 
intermediate phase characterized by angle 9, per unit cell, for the case of classical spins is obtained as: 


£7(0) = -2 S i S 2 (F cos 9 + K sin 0) - A [s \ + 5 1) (2) 

in terms of the magnetic coupling parameters F, K, and A. 

Figure[3ji shows a phase diagram in (F, K) space, obtained by comparing the classical energies of three phases: 
the collinear (0 = 0), the orthogonal (9 = n/2) and the 9 — 7t/ 4 intermediate phases for A — 1 and equal spin 
magnitudes S \ -Si- 

To explore the dependence on angle 0, one can also compare analytically the energies given by Eq.(0 for a 
generic 0-angle intermediate phase, with the energies of the collinear and the orthogonal phases. The two bound¬ 
aries for the three observed regions in (F, K) space (like shown in Figure[3ji, for 0 = 7 t/ 4) are found to be given by 
linear functions, with respective 0-dependent slopes as shown in Figure [3]) : 


m i 


1 - cos 0 
sin0 


m 2 


cos 0 
1 - sin 0 


(3) 


independent of the sublattice spin magnitudes. Thus, fixing the angles of the phases included, the phase diagram in 
(F, K) space would not be modified even if different spin values were used, as the energies of the different phases 
are rescaled proportionally. In Figunj3}i the angular dependence of m \ and mi is exhibited. 

From ([2} one can also obtain the angle 0 mm which leads to the intermediate phase with minimum classical 
energy, for any set of coupling parameters: 


dE[(9) 

d9 


0 m i n = arctan 



(4) 
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Figure 3: Analysis of stability of the different phases in the classical model. Parameters: A=1;5 i= 52 = 0.5. 
(a) Classical phase diagram including: the 6 = n/4 intermediate phase, the collinear phase (9 = 0) and the 
orthogonal (6 = n/2) phase, (b) Angular dependence of the slopes of phase boundaries m\ and m 2 given by Eq.©- 
(c) #min as a function of F, for K = 0.3 (dashed line); 0 m ; n as a function of K, for F = 0.3 (dotted line). 


Notice that <9 m i n is independent of NNN coupling parameter A, and only depends on the “intra-dimer coupling 
ratio: K/F. In Figure [3]: we plot /9 lmn along specific lines in parameter space, marked in Figure [3}i. We show the 
monotonously decreasing 8 mm as a function of F, for K = 0.3 and, as could be expected, confirm that for F = 0 the 
orthogonal phase 6 = 7r/2 represents the stable ground state. The monotonous increase of 0 m j n with K, for F = 0.3, 
is also shown, and we here confirm that for K = 0 the collinear phase 8 = 0 is stable. Notice also, in Figure[3}:, that 
the 8 = 7t/ 4 intermediate phase will only be stable when K = F. 

Finally, in Figs. [4a] and [4b] the classical energies of the collinear phase, the orthogonal phase and the 8 = /9 lmn 
intermediate phase, are plotted as functions of F and K. It becomes clear that, with our simplified model, for each 
(F, K) set of parameters, one intermediate phase (the one with 8 - f> ulm ) is always the classically stable ground 
state, merging with the proposed collinear and orthogonal phases in the appropriate limits. 



(a) Dependence on F, at K = 0.5. 


(b) Dependence on K, at F = 0.5. 


Figure 4: Classical energy of the: collinear, orthogonal and 8 = 8 m i„ intermediate phases as a function of coupling 
parameters F in (a), and K in (b) Other parameters: A = 1;5'i=5 , 2 = 0.5. 
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3.2. Calculation of the quantum spin excitations 


To calculate the quantum magnons of our chain model for nickelates at low temperatures, given by Hamiltonian ©, 
we start by performing a local rotation of the S 2 sublattice spin quantization axes by an angle 8, with respect to the 
S 1 sublattice, using the following transformation (S y n - S' n ): 

S* = cos6S* -sin#S~ S z n - sin OS* + cos0S] (5) 

Next, with the Holstein-Primakoff transformation the Hamiltonian is rewritten in terms of bosonic operators, 
and the Linear Spin Wave approximation (LSW) is used: 


I S n , 

Vt ( o ' 

(6) 

\Sn ,, f . , „V ■ jSn ,pf , N S n -—S n +b n b n 
sj — (Pn + bn) S„ = -I y — (b' n - b„) 

where (a\ a) operators refer to the spin up sublattice, and (/;’, b) to the spin down sublattice. Though the resulting 
Hamiltonian includes one, two and three operator terms, in LSW we consider only those written in terms of two 
operators. One operator terms change the energy of ground state. Introducing the Fourier transform of the boson 
operators, we obtain the following Hamiltonian: 




"H =} {A\ (a^ a, +b'' b. ) + /L (fl ^ a, +b' 1 b, ) + <5+ (fl 1 a, + a, a^ + Z? 1 b, + b. /? 1 ) 

Z_i [ 1 \ 1,9 >'9 3,9 1 v 2., 2.9 4,5 *.1J + \ 1, 9 2,? 1.9 2., 3 ,, 4,9 3 ,, i<q J 

1 

+ { a t a L + fl >., a 2 .-, + bl q bl q + b, q b 4 _ q ) + ffq)a}bl + ri {q)a lq b 3q 

+ yliq) a t h l + yiiq) a^b^ 


(7) 


with the coefficients defined as (s = PS 1 S 2 )■ 

A 12 — 5 2,1 cos $ + K sin 6j + 2A S b± - -^-JF(cos 8 + 1) + K sin jipiq) = A 5 1,2 (l + e ,?a ) (8) 

Finally, by paraunitary diagonalization of the Hamilonian, we determine the energies of the magnon excitations. 


3.2.1. Magnon predictions for the different phases 

In the following, we discuss the quantum magnons obtained for Hamiltonian ([7} using different pairs of (F, K) 
couplings, as detailed in figure|5a] The AFM coupling between “dimers”, A = 1, was considered as unit of energy, 
and except for Fig.[5b]the spin magnitudes where chosen as S \ - S 2 - 0.5, which corresponds to 6 = 0. 

In all cases we obtained 16 eigenbands, corresponding to the 16x16 Hamiltonian matrix, of which only the 8 
positive branches describe magnon excitations. 

First, notice that with the different pairs of (F, K) couplings used, in Figs. l5bll5fl we are exhibiting the quantum 
magnons which correspond to the intermediate phase with lowest classical energy corresponding to angles /9 lmn 
ranging from 0 to tt/2. 

Now, notice that the lowest magnon branch both in Figures [5b] and [5f] is identical to the spin excitation of an 
antiferromagnetic chain. In Fig. [5f| the lowest branch corresponds to the collinear phase (6 = 0) with parameters: 
F = 0, K = 0.5, hence in © the effect of K on the magnons disappears, due to the sin 8 factor. Thus, the only 
non-vanishing coupling term in the Hamiltonian is determined by A, hence the antiferromagnetic chain excitation 
obtained. In contrast, in Fig. [5b]the lowest branch corresponds to the orthogonal phase (8 = n/2) with parameters: 
F = 0.5, K = 0. In this case, in Eqs. © the effect of F does not disappear as K in the previous one. Nevertheless 
we checked that the dispersion of an antiferromagnetic chain results even when F —» 00 . This can be understood 
because in this limit the system behaves as an antiferromagnetic chain composed by dimers, each of them consisting 
of two consecutive spins coupled by F. 
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(a) Points of (F, K) space, where the predicted 
quantum magnons are shown next. Colour-depth 
indicates the value of classical angle 0 m j n , as in¬ 
dicated by the scale included. Unless otherwise 
stated: A = 1, $ i = S 2 = 0.5 



(c) Magnons corresponding to: F = 0.5, K = 0.3 (A in 
Fig|5al. Here, the intermediate phase that minimizes the 
classical energy corresponds to angle: 0 m i n = 30.9°. 



(e) Magnons corresponding to: F = 0.3, K = 0.5 (♦ in 
Figl5al. Here: 0 min = 59°. 



(b) Magnons obtained with intra-dimer couplings: 
F = 0.5, K = 0 (■ in Fig.[5a}, for different phases. Solid 
line: collinear phase (0 = 0 = 0 nllI1 ). Dashed line: or¬ 
thogonal phase (8 = n/2). Dotted lines: collinear phase 
(0 = 0 = 0 m i n ) for S i = 0.6 and S 2 =0.4. 



(d) Magnons corresponding to: F = 0.5. K = 0.5 (★ in 
Fig|5aJ. Here: 0 min = 45°. 



(f) Magnons corresponding to: F = 0, K = 0.5 (• in 
Fig|5al. Orthogonal phase (0 = n/2 = 0 m i n ): dashed 
lines: for S 1 -Si - 0.5; dotted lines: for ,S 1 = 0.6 
and S 2 = 0.4. 


Figure 5: Magnons for RMO 3 nickelates obtained with our spin chain model for different parameters. Inter-dimer 
AFM coupling A = 1; spins: S \ - S 2 — 0.5, other parameters as detailed in subcaptions. Unless otherwise stated: 
collinear phase with 0 = 0 (solid line); orthogonal phase with 6 = n/2 (dashed lines); 0 m in-intermediate phase 
(dot-dashed lines). 
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As shown in Figure|5b] the magnons obtained for the collinear phase (solid line), which is the classically stable 
phase, are all degenerate. We checked that these excitations correspond to spin-flips between two NNN inter-dimer 
spins: i.e. in sites 1 and 3, or 2 and 4 with the notation of Fig. [2b] Both of these spin flips involve the same energy 
cost related to coupling A, thus being degenerate. As one would expect, we find that considering a certain degree 
of charge disproportionation, i.e. S \ ^ S 2 in the model, some degeneracies are broken due to the lower simmetry 
of the system: we exemplify this in figure [5b] by including the case S 1 = 0.6, Si - 0.4 (6 - 0.2), plotted with 
dot-dashed lines. This effect is largest at the Brillouin zone edge X. 

Regarding the orthogonal phase. Figure [5fl shows its excitations for F = 0 and K = 0.5 (• in Fig. [5aJ), case in 
which it corresponds to the classical stable phase. Notice that turning on the K coupling also reduces the symmetry 
of the system and, as previously mentioned, breaks magnon degeneracies. Though the analysis of the obtained 
excitation modes is more complex for this phase, the symmetry breaking of the two highest magnon branches can 
be understood as follows. The highest energy magnon branch (with energy ~ 0.5 at T, in Fig. [5]} involves spin-flip 
excitations between two NN inter-dimer spins in sites 1 and 4, and in sites 2 and 3, involving an energy cost related 
to A and K couplings. A lower energy cost is payed exciting the magnon branch below it (with energy ~ 0.4 at F, 
in Fig. [5TJ>: which we checked corresponds to spin excitations of two NNN inter-dimer spins (i.e. between spins 1 
and 3, and spins 2 and 4), in which case only the K coupling is affected. 

It is also interesting to compare the effect of the charge disproportionation in the two cases exhibited in Figs. [5b] 
and[5TI respectively. While, as discussed above, in the collinear phase of Fig. [5bl degeneracies are clearly split by 
6, we observe that in the orthogonal phase depicted in Fig. [5f| no new degeneracy splittings appear, in addition 
to those originated by the presence of the DM coupling K. In fact, the main ^-dependent effect produced by 6 
in the orthogonal phase is to increase the size of the magnon gap between the upper and lower pairs of branches 
originated by K. According to our results, the different numbers of magnon branches observed might thus be used 
to distinguish between the collinear and orthogonal phases, and even to quantify the charge disproportionation. 

Also, the figures [5b][53 show that in the parameter ranges considered none of these phases becomes unstable, 
being all excitation energies positive. Even though this does not allow us to indirectly determine which of the 
studied phases would represent the ground state of our chain, it is plausible to infer that the quantum ground 
state should be similar to the intermediate phase that minimizes the classical energy (i.e. / flmin : plotted with dot- 
dashed lines in Figure 0). This might be justified observing that, in figures l5b[|5TI the 9^ -intermediate phase 
has excitations with higher energy than the other phases, thus it seems more difficult to create excitations and 
eventually destabilise the 0 m ; n -intermediate phase. 


4. Summary 


A simplified localized spin chain model was proposed to study the generic intermediate phase in nickelates, 
able to describe a variable charge disproportionation and relative orientation of consecutive spins along the chain. 

The model includes the following magnetic couplings: nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
Heisenberg-like interactions, respectively for the ferromagnetic and antiferromagnetic couplings present in the 
collinear phases. To describe the non-collinear phases, we also consider a NN Dzvaloshinskii-Moriva lEllITTtl -tvDe 
coupling to allow for the possibility of a relative angle between the different magnetic sublattices. 

We studied: (i) the classical stability of the collinear, orthogonal, and intermediate phases, as possible ground 
states for these compounds; and: (ii) the quantum ground state indirectly, by calculating the spin excitations 
resulting from each of those phases. 

Our classical results show that for each set of NN (intra-dimer) ferromagnetic and DM magnetic couplings, 
always an intermediate phase characterized by an angle 0 m j n , corresponds to the most stable classical ground state. 
From the measurements of the two Ni spin magnitudes by Fernandez et al. 112j], in terms of our model one can 
obtain the following estimation for the relative orientation of consecutive spins along the chain, 6 ~ 80 c "' c , and 
intra-dimer coupling ratio: K/F ~ 6. 

Regarding the quantum mag netic excitations: with our simplified model, we predict the spin excitations to 
be expected for the collinear QT12] and the orthogonal phases Tit 18 ] so far proposed for these compounds, as 


well as those of the intermediate phase. For the collinear and orthogonal (6 — tt/2) phases, we predict differences 
in the magnon spectrum which would allow to distinguish between them in future inelastic neutron scattering 
experiments. In particular, the number of magnon branches would differ between these phases, and the charge 









disproportionation present might also be quantified: either by the number of branches in the collinear phase, or by 
the size of the magnon gap in the orthogonal phase. 

Our present study represents a first step towards an understanding of the complex three-dimensional ground 
state of rare-earth nickelates, to enable comparison with the results of future inelastic scattering experiments. The 
latter are especially desirable in these compounds, since the available neutron diffraction data could not discrimi¬ 
nate between the different ground states proposed. Material specific ingredients, as well as the different possibilities 
of three-dimensional stacking should be considered in future research work. Also, if no obvious signs of destabi¬ 
lization of any of these phases are observed in the predicted magnons, a direct study of the quantum phase diagram 
would be important. 
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